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Abstract 

The proton-proton fusion process pp — * de + v e is calculated at threshold to 
fifth order in pionless effective field theory. There are two unknown two- 
body currents contributing at the second and fourth orders. Combined with 
the previous results for u e d and V e d scattering, computed to third order in 
the same approach, we conclude that a ~10% measurement of reactor v e d 
scattering measurement could constrain the pp — > de + v e rate to ~7% while a 
~3% measurement of u e d — ► e~pp could constrain the pp rate to ~2%. 
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The reaction pp — > de + v e is of central importance to stellar physics and neutrino astro- 
physics. Bethe and Critchfield proposed it to be the first reaction to ignite the pp chain 
nuclear reactions that provided the principal energy and neutrinos in the Sun P|. Major 
efforts have been made to provide precise theoretical predictions for this reaction at zero en- 
ergy |2|-|T3|. No direct experimental constraint is available for this process, and the accuracy 



of the theoretical predictions must always be weighed against the availability of empirical 
constraints. However, Schiavilla et al. recently calibrated the matrix element with comple- 
mentary calculations of tritium beta decay to obtain an estimated uncertainty of less than 
1% in the potential model. 

Alternatively, effective field theory (EFT) can provide a connection between pp fusion 
and other processes that might be accessible experimentally. Specifically, there is a direct 
connection between the reactions pp — *> de + v e and v e d — > e~pp in that they both involve 
the same matrix elements. Further, in EFT, they can both be shown to depend on a 
single unknown counterterm (or two-body current) at third order. Otherwise, the nuclear 
physics input is identical. Further, all four v(y)-d breakup channels depend on this same 
counterterm, meaning that a measurement in any one channel implies a measurement in all 
channels - including pp fusion. To date, measurements are available for V e d breakup using 
reactor antineutrinos with ~ 10-20% uncertainty [p-JTS] and possibly ~ 5% in the future 



Also, the ORLaND proposal could provide a measurement of v e d — ► e pp to a few 



percent level []19 |. It is important to understand the constraints that these measurements 
could on pp — > de + i , e . 

At the present time, v e d and V e d breakup processes have been studied to the third order 
in pionless EFT depending on the unknown two-body counterterm L\a- Through varying 



Li : a, the four channels of potential model results of refs. [£0J and |H] are reproduced to 
high accuracy. This confirms that the ~ 5% difference between refs. [2D] and |2TJ is largely 
due to different assumptions made to short distance physics. We will discuss this difference 
later in this letter. 

For pp fusion, an pionless EFT calculation has been performed to second order by Kong 
and Ravndall p2| , with a resulting dependence on L\ t A- In this letter, we use the same 



approach to push the calculation to fifth order, which introduced another unknown two-body 
counterterm (Ki } a) that first enters at fourth order. Constraining K^a using dimensional 
analysis, we conclude that a measurement of reactor V e d with 10% uncertainty can constrain 
pp fusion to 7% while a measurement of v e d — > e~pp to 3% would constrain pp fusion to 2%. 

The method we use is pionless nuclear effective field theory, EFT (ft) P5| , treating the 
electromagnetic interaction between protons in the manner developed in @,[24[] . The dynam- 
ical degrees of freedom are nucleons and non-hadronic external currents. Massive hadronic 
excitations such as pions and the delta are integrated out, resulting in contact interactions 
between nucleons. The nucleons are non-relativistic but with relativistic corrections built 
in systematically. Nucleon-nucleon interactions are calculated perturbatively with the small 
expansion parameter 

_ (1/dpp, qf,p,aM N ) 

A [) 

which is the ratio of the light to heavy scales. The light scales include the inverse S-wave 
nucleon-nucleon scattering length l/a pp (a pp = —7.82 fm) in the 1 S , , pp channel, the deuteron 
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binding momentum 7(= 45.69 MeV) in the 3 Si channel, the proton momentum in the center- 
of-mass frame p, and the fine structure constant a(= 1/137) times the nucleon mass M^. 
The heavy scale A, which dictates the scales of effective range and shape parameters is set 
by the pion mass m n . The Kaplan-Savage- Wise renormalization scheme |25| is used to make 
the power counting in Q transparent. This formalism has been successfully applied to 

many processes involving the deuteron [23 , 27| including Compton scattering p8|j29|| , np — ► 
dj for big-bang nucleosynthesis vd scattering for physics at the Sudbury Neutrino 

Observatory [Q, and parity violation observables |34| . 



There are other power counting schemes which yield different orderings in the pertur- 
bative series and each has certain advantages. For example, the z-parametrization |35| 
recovers the exact deuteron wave function renormalization at second order; the dibaryon 



pionless EFT |28] resums the effective range parameter contributions at first order and sim- 
plifies the calculation tremendously by cleverly employing the equations of motion to remove 
redundancies in the theory. While one of these power countings could lead to more rapid 
convergence in any given calculation, for a high-order calculation as we present here the 
distinctions between different expansions are negligible. 

Ignoring for the moment the weak interaction component, the relevant lagrangian in 
EFT(Tf) can be written as a derivative expansion 



C = rf id + 



2M 



N 



N , 



i 



Si 



c 
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(N T PiNy(N T V 2 PiN) + h.c. 
- (N T V 2 P i N)\N T W 2 P l N) 



(N T V 4 PiNy(N T PiN) + h.c. 



32 

3 Si — > 1 S'o, P{ — > PA + 



(2) 



where V = ^? — V and where Pi = o^o^/v^ and Pi = <J2 T 2 r i/V8 project out 3 Si and 
1 S'o channels respectively with er(r) acting on spin(isospin) indices. The coupling constants 
have been fit to data in p3|[32|. We only perform the calculation at threshold and thus 



the relativistic corrections and momentum transfer (q) effects are suppressed by factors 

Vo. Thus, with the goal of presenting a 
the nonrelativistic and zero recoil limits are 



of 7 2 /M|r and q 2 /7 2 , both of which are < 1%. Thus, 



calculation with a precision of less than 1, 
suitable approximations for us to make. 

The weak interaction terms in the lagrangian can be written as 



C 



w = _ ^§ Z + J M +h - c - + 



(3) 



where Gp is the Fermi decay constant and 1+ 
hadronic current has vector and axial vector parts 



z^7 M (l — 75 )e is the leptonic current. The 
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The vector current matrix element vanishes in the zero recoil limit since it yields terms 
proportional to the wave function overlapping between two orthogonal states. The time 
component of the axial current gives center of mass motion corrections which vanish in our 
approximation. The spatial component of the axial vector is the dominant contribution, and 
can be written in terms of one-body and two-body currents 

K = 9 ^r-a k N 



2 

+L 1>A 

K ltA 
8 

+ ••• , 



(N T P k N)\N T P~N) +h.c. 
(^N T V 2 P k N^ (iV T P~iV) + (N T P k Ny (iV T V 2 P~iV) + h.c. 



(4) 



where qa = 1.26, t~ = {t\ — ir 2 ), and the values of L^a and K\ a are yet to be determined 
by data. 

The hadronic matrix element is usually parametrized in the form 



d;j 



pp 



= 5uCV 



32vr 



where j is the deuteron polarization state, and 



Cn 



2nr] 



aM, 



N 



e 27 "> - 1 ' 



2p 



(5) 



(6) 



is the well-known Sommerfeld factor. In the center of the Sun, the scale of p is ~1 MeV. 
Thus the Sommerfeld factor C 2 changes rapidly with respect to p while A(p) does not. 
It is sufficient to keep the p 2 correction for A(p), since the higher order correction is 0( 



(p/aM^ 



N 



Using e to keep track of the Q expansion we find that, to fifth order in the Q(e) 



expansion, A(0) can be written in the compact form 
1 



A(0) 



e x - 7a pp [1 - xe x E 1 (x)] - e-f a 



pp 



e 2 7 2 _ 
L\a — KiA 



l + O (e 5 ) . (7) 



This expression can be expanded to e 4 , and then e should be set to 1. p d = 1.764 fm is the 
effective range parameter in the 3 Si channel, x — cnM N /j and 



roc g 

E 1 ( X ) = / dt- 

Jy I 



(8) 



Li a and K\ a are the renormalization scale /^-independent combinations of the /x-dependent 
parameters L x a, K\ a and the nucleon-nucleon contact terms C 2 and C 4 



Li,A 
Kl,A 



- (/i - 7) 


Li t A 




. 9 a 


- {n - 7) 


K\,a 




. 9 a 



M 



N 



(c^ + a 



(d) 
2,-2 



— — M N (C^\ + 2Cf 3 ) 



(9) 
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^From ref. |3~2" |, we have 

r ( PP ) i!L I 

^ p -^ + N { ^ + ^ E ) 

n {p P ) _ M N ( P p) r ,(pp)2 Mpp) _ Mn (pp) r ,(pp)2 

°2,-2 — g^T'O °0,-l ) °4,-2 _ 47T 0,-1 1 

{d) 2tt p d {d) _ 7T p 2 d 

C2 - 2 "^(^f' C4 '- 3 " M*(^?' (0) 

where the second subscripts of the C's denote the scaling in powers of Q, 7^ = 0.577 is 
Euler's constant, and = 2.79 fm is the pp channel effective range. We take the pp 
channel shape parameter to be the same as that in the np channel = —0.48 fm 3 . Any 
errors introduced by this are small numerically. 

After expanding to e 4 and setting p = m n , we obtain 

A(0) = 2.58 + 0.011 (^) - 0.0003 (^) . (11) 



At p = m^, dimensional analysis as developed in refs. p3| , p5| would favour 

1 „ 

\ l i,a\ ~ ; Z2 ~ 6 fm > 

m n (m n - 7) 

\K 1A \ « ^ « 20 fm 5 . (12) 

m£ (m^ - 7) 

Two observations follow. First, if we take these naively estimated values (with positive signs, 
for example), then the expansion of A(0) converges rapidly 

A(0) = 2.51 (1 + 0.039 + 0.029 - 0.010 - 0.0001) . (13) 

For other sign combinations, the series also converges rapidly and higher order effects are 
<C 1%. Second, eqs. (PD and ( |T2"D show that K\^a is likely to contribute to A(0) at a level 
less than 1%. Thus eq. ([TT|) is precise to 1% even with K 1A set to zero, meaning that we 
can write 

A(0) = 2.58 + 0.011 (j^g) + 0(1%) . (14) 

This the central result of this paper. 

Ultimately the value of L\ a must be extracted from experimental data. Before we 
address this issue, let us look at what values of L 1A are found in fits to potential model 
calculations. Using the results of ref. [ 32 1 , the recent potential model results for u(v) — d 
breakup of Nakamura, Sato, Gudkov, and Kubodera (NSGK) []2T] are equivalent to 

L™ GK = 5.6 ±2fm 3 , (15) 

while the results of Ying, Haxton and Henley (YHH) |2D| are equivalent to 
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Lia_ H = 0.94 ± 2 fm 3 . (16) 



The uncertainties represent a conservative estimate from EFT of 3% from ref. (even 
though the NSGK results can be reproduced within 1% using the central value of £i,a)- 
Based on the size of the third-order contribution, the actual uncertainty may be as small 
as 1%, but further analysis is required to ascertain this. The (undefmable) errors from 
the potential models themselves are not included. As mentioned before, the differences 
between these two calculations are in their treatment of two-body physics. NSGK uses 
a model to include axial two-body meson exchange currents, while YHH includes vector 
two-body currents but not the more-important axial two-body currents. Given that Lx,a is 
representative of the dominant axial effects in EFT, it is not surprising to see substantial 
differences in the value inferred by each calculation. 

Similarly, eq. ([14]) translates the pp fusion rate (A 2 (0) = 7.05-7.06) calculated by Schi- 
avilla et al. into 

L Schiavilla et al. = q 5 ± 3 4 fm 3 ? 

which is consistent with NSGK. 

Experimentally, it is easy to relate the EFT v(v)-d scattering results and pp fusion rate 
through the third-order results for v{y)-d provided in [j32]| . The results for each channel are 
parameterized in the form 

a{E v ) = a{E u ) + b(E u )L 1<A + 0{< 3%) , (18) 



where a and b are functions of neutrino energy. These results, tabulated in ref. [j52"|, can be 
easily related to the pp fusion rate through eq. QT4]). 



For reactor V-d scattering, one expects the rate to be peaked around 8 MeV, which can 
be interpreted to mean a measurement precise to 10% can determine the pp fusion rate 
(oc A 2 ) to 7%. With the experimental precision further improved to 5%, the pp fusion rate 
can be determined to 4%. Alternatively, the proposed ORLaND |19[ detector might measure 
v e d — > e~pp to 3%, in turn constraining the pp fusion rate to 2%. 

In strong interaction processes involving external currents, delicate relations between op- 
erators are required to guarantee that there are no off-shell ambiguities in the final results. 
For example, in two-body systems two-body currents serve to absorb the off-shell effects 
from the two-body strong interactions. This means that different models that reproduce the 
same on-shell nucleon-nucleon scattering data might need quite distinct two-body currents 
to deal with off-shell effects. In the case of pp fusion or v-d scattering, electromagnetic 
matrix elements cannot be used to constrain weak matrix elements. Their operator struc- 
tures are quite different at the quark level even they might appear the same in hadronic 



level. Using tritium /3-decay to constrain the two-body current |L2| is an excellent idea, in 
principle. However, contributions from three-body currents that are required to absorbed 
three-body off-shell effects are not yet constrained. In light of these facts, the need for a 
precise experimental measurement of u e {p e )-d breakup cannot be overstated. 
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